A coherent and passive one dimensional quantum memory 
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We show that the state of a flying qubit may be transferred to a chain of identical, (near) ferromagnetically 
polarised, but non-interacting, static spin- 1 particles in a passive way. During this process the flying qubit is 
coherently polarised, emerging in the direction of the majority static spins. We also show that this process 
is reversible for at least two flying qubits injected sequentially and thus has the potential to be exploited as 
a passive quantum memory to encode the flying qubits without the necessity of resetting between successive 
encoding operations. We show that the quantum information may be spread over many static spins in the 
memory chain, making the mechanism resistent to spin decoherence and other imperfections. Among some 
potential architectures, we discuss implementing the memory in a photonic waveguide embedded with quantum 
dots, which is resilient to various possible errors. 



Robust quantum state transfer (QST) plays an important 
role in the field of quantum information processing (QIP), 
achieving quantum transmission of data through space or 
time Over the past decade numerous efforts have been 
made in this area, and many potentially feasible approaches 
have been suggested for both state transportation 12[] and stor- 
age Hi in a variety of physical systems. However, in most 
instances, exercising the required level of active quantum con- 
trol remains a challenging aspect of the current technology: 
errors are most likely introduced. On the other hand, it has 
recently been proposed that schemes utilising iterative appli- 
cations of quantum maps can perform certain QIP tasks with 
reduced level of quantum control @] ■ 

In this Letter, we work with an iterative setting for a coher- 
ent quantum memory under limited quantum control; more 
specifically, we propose a scheme to sequentially transfer the 
states between a number of flying qubits and a long (memory) 
chain of TV identical, ferromagnetically polarised, but non- 
interacting, static spins in a passive way (see Fig.UJ. 
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FIG. 1. Illustration of the passive memory system (not to scale) with 
the associated QST scheme. The coherent quantum memory consists 
of a sufficiently long ferromagnetic chain \F) of non-interacting, 
static spins (red). Each flying qubit (blue) enters the front end of 
the chain, interacts with the static spins sequentially and eventually 
emerges as polarised t)/ at the tail end. The chain can then en- 
code subsequent flying qubits sequentially in the same fashion. To 
read out from the memory, one simply injects a polarised qubit |t)y 
back from the tail, and the state of the last encoded flying qubit is 
then recovered automatically out from the front. More states can be 
recovered sequentially in the same way by further injecting back po- 
larised qubits |t), from the tail. There is no need to reset between 
successive encoding (decoding) rounds. Any multi-partite entangle- 
ment between injected qubits is also recovered in the read operations. 



We start by considering the first flying qubit in some ar- 
bitrary state l^i)^ = cx-x |t) + Pi It) an d its associated 
QST process with the static spins in the ferromagnetic chain 



l^) c = It. 



As shown in Fig. [TJ during encoding 



(or the 'write' operation) \i/>i) * ± enters from the front into the 
chain and interacts with the static spins si to % sequentially. 
We model the interactions by an effective Hamiltonian cou- 
pling each flying qubit and the k th static spin of the following 
form: 



— (o- f o- Sk +a f a Sk ) 



g k (o- f + a s _ k +a J _o- s + k ), (1) 



where o± = (a x ± ia y )/2 are Pauli spin-flip opera- 
tors. The g k are the XY exchange coupling strengths 
that depend on the separation of the qubit and the spin 
s k , and are hence time dependent due to the mobile qubit. 
Thus for a general state |# fc (i)) = U k (t) \$> k {0)) of 
the qubit and the spin s k , the time evolution operator 



U k (t) 
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) , where 6 k (t) 



Jo 9k(t')dt' /h is constant when [0, t] is chosen so that g k (0) 
and g k (t) are negligible. In the basis |t/ts fc )> lt/tsj> 

|4-/t«k)> |t/ts fc ), 
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For the Heisenberg model with an extra term ^-a{a s z 
Eq.[TJ U k takes the same form apart from an extra phase factor 
of e l9k for each trigonometric term. Let us first consider the 
XY case with 9 k = 9 e (0, tt/2] V k. 

Starting from the total state \ipi,F) = \ipi) ^ \F) C , we ap- 
ply U k as in Eq. |2] to the flying qubit and the k th static spin 
sequentially for si through to sn- The state |t, F) with am- 
plitude ai remains the same during this write operation, while 
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||, F), with amplitude fti, evolves as follows: 

\i,F)=>i cos9 ||, F) - i sin |t,F) 

^ 2 cos ( cos 6» ||, F) - i sin 5j |f, F)) - i sin ||, F) 
=^ fc ...Vfc = 3,...,iV-l 

n <zs iV — y oo 

^jv -isin0^cos fc - 1 5^ ||,F) +pe^0||,F) (3) 
fc=i 

where =>fe corresponds to the occurrence of an interaction 
event for the flying qubit with the static spin su, and de- 
notes a spin down-flip in the k th position of the chain. Note 
that the only component which evolves further in each super- 
position is underlined in Eq. [3] Combining both parts, we see 
that the flying qubit emerges as polarised \t) f 1 with probabil- 
ity 1 — cos 2N — > 1 for N —> oo. In this limit, the initial 
quantum information which \ipi) ^ held before encoding has 
now been transferred to the chain whose collective state reads 
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This is model independent by conservation of total spin, and 
so the argument is equally valid for the Heisenberg model (up 
to factors of e lke ). Note that the probability of a single down- 
flip at site k is |/?i | 2 sin 2 cos 2 ( fe-1 ' 0, a quantity that decays 
exponentially along the chain. Summing over all these proba- 
bilities gives | /3i | 2 , as expected by conservation of total prob- 
ability. 

To read out the state \ipi) of the original flying qubit from 
the memory at a later stage, we inject a polarised flying qubit 
|t) . from the tail back to the chain, i.e., in the opposite di- 
rection as for the encoding operation (see Fig. [TJ. This is the 
simplest decoding method for the memory system, especially 
when the level of quantum control is limited. The flying qubit 
sequentially interacts with the static spins through % to si 
and, provided it has the same kinetic energy as the emitted 
qubit /i, the total state of the system after decoding becomes 

OO \ 

at |t, F) - /3i sin 2 ^ cos 2 ^" 1 ' ||, F) J - ifr sin x 
fe=i ' 




\t,F) .(5) 



This is straightforwardly proved by applying the two-spin uni- 
tary operations sequentially as detailed in the Supplementary 
Material ifTTll . 

Thus, after the read operation the chain returns to the orig- 
inal ferromagnetic state and is disentangled from the flying 
qubit, which now emerges from the front as 



In other words we have recovered the original state of this fly- 
ing qubit, up to a phase flip which can be corrected by a simple 
er| gate. This will always be the case for the XY model, and 
importantly distinguishes the decoding process from a time 
reversal operation. However, the same read operation does 
not work for a Heisenberg-type coupling (except for special 
cases), since the extra phases e lk 9 invalidate the cancellation 
inEq.0 

To see how the size N of the memory chain required to 
store one qubit scales with the coupling strength 0, we note 
that the only condition is cos 2W — s- for sufficiently large 
N. Given any error tolerance e > such that cos 2N < e, we 
require 



2N(6) > 
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Thus we see that the required memory size increases signif- 
icantly with decreasing coupling strength for fixed e. For 
example, with e ~ 1CP 4 and = 1, iV m j n ~ 8 rising to 
~ 920 for = 0.1. However, in the weak coupling regime, 
the qubit state stored in the chain is delocalised; each static 
spin only shares a fraction of the total information, and thus 
has the potential to store more. On the other hand, for strong 
couplings a relatively short chain can already store a num- 
ber of flying qubits. In the special case where = 
one polarised static spin in the chain is sufficient to store 
a flying qubit, since now Uk in Eq. [2] is simply a SWAPk 
gate (up to a phase of '— i'). Encoding of subsequent fly- 
ing qubits is thus also obvious, as are their readout proce- 
dures (since SWAP 2 , = I up to a phase of '-1'). These qubits 
are stored locally in the memory: By induction and linear- 
ity, one finds that n ordered flying qubits \ip n )f ■■■H'i) f 1 are 
encoded sequentially into the memory and the chain state be- 
comes \ip n ) Sl ■■■ \ipi) s It ■■■ t) s after the whole write oper- 
ation. Once read back, we obtain \ipn)f ■■•\' t Pi)f 1 > where 
~ denote the appropriate phase changes for the XY model. 
For this special case the argument also holds for the Heisen- 
berg exchange model, though the extra phase flips will not 
be present. Indeed in the Heisenberg case the decoding pro- 
cess is nothing but a time reversal of the encoding operation. 
In fact, realisation of time reversal (0 — > —0) for arbitrary 
would make the memory system viable even for the Heisen- 
berg case in general. However, this would require tuning from 
ferromagnetic coupling gu (write) to anti-ferromagnetic — gh 
(read), which is impractical. 

We now consider the effect of decoherence on a chain in 
which one qubit state |+) * = ( \t) f t + |4-)/j ) /V2 is stored. 
Both the write and the read operations are assumed to be fast; 
the principal effect of decoherence is on the static spins for 
the time (t) during which the qubit is stored in the chain. We 
model this using the Lindblad master equation tZD, and by re- 
stricting to the subspace consisting of only zero or one excita- 
tions among the static spins. We simulated the behaviour for 
the memory chain of a fixed length TV = 100 under dephasing 
errors. We also vary such that the qubit is stored in the first 
N s (< N) number of static spins with e ~ 10~ 2 (see Eq.|7]), 
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to monitor the effects of spreading the quantum information 
stored in the memory. For homogeneous dephasing, in which 
each spin is subject to the same, independent decoherence pro- 
cess, the total decoherence for the memory is essentially the 
same for all 9 values. In our simulations we take a dephasing 
rate T = 1 MHz for each static spin and each resulting curve 
coincides exactly with the red curve in Fig. |2ja), regardless 
of the chosen N s < N. In other words, the decoherence rate 
does not depend on how local or distributed the information is 
in the quantum memory, since the relevant qubits which con- 
tain the quantum information each decohere at the same rate 
r. It also follows that the fidelity must eventually saturate to 
i when all the memory qubits lose their quantum information. 



Fidelity 
1.0 
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FIG. 2. Plots for the fidelity of the retrieved qubit, relative to the input 
1+)^, against the storage time r under inhomogeneous dephasing 
with rate F — 1 MHz (coherence time 1 fj,s) for the fc* static spin 
and zero for all others (iV = 100). We vary 9 such that the qubit 
is stored in the first N s number of static spins with e ~ 10~ 2 (see 
Eq.[7]l: (a) The qubit is stored entirely in the first qubit (9 — tt/2); (b) 
The whole chain of 100 spins store the qubit collectively (9 ~ 0.30). 



On the other hand, for inhomogeneous decoherence pro- 
cesses in which the spins decohere at different rates, distribut- 
ing the quantum information does reduce the variance of the 
total memory decoherence rate compared with locally stored 
information, and this difference can be very large. We illus- 
trate this in Fig. |2] for the case when just one static spin is 
subject to decoherence (for example due to the proximity of a 
magnetic impurity). When the quantum information is stored 
locally (in the first qubit with 9 = ^ in Fig. [2ja)) then the 
fidelity F saturates to i when the decohering spin is the first, 
whereas F = 1 if the decohering spin is any of the other N— 
1 spins. However, when the quantum information is spread 
over all N spins (9 ~ 0.30 in Fig.|2jb)), the saturated fidelity 
becomes > 1 — ^ , independent of which spin is subject to 
the decoherence. This result may easily be extended to cases 
when more than one spin decoheres, or the spins decohere at 
different rates. Spreading the information over many qubits 
will again smooth the statistical fluctuations in information 
loss. 

The general result for encoding and decoding an arbitrary 
number of input flying qubits with arbitrary multi-particle 
superpositions remains presently a conjecture, but we have 
strong evidence that it is true. In particular, we have extended 
the proof from one to two flying qubits, including entangled 
states, as well as for the case where there is only one |4)* 



among any number of flying |t) * qubits. Both of these proofs 
can be found in the Supplementary Material | l lfl, but we out- 
line some essential ingredients of them here. 

For the second of the two proofs just mentioned, we define 
a th collective 1-spin down-flip operator on the chain \F) 
arising from Eq. [3] 

oo oo 

= $k ■■= EH sin0cosfc "^) § fc • (8) 



fc=l 



k=l 



We find J2k \ a< o\k)\ 2 = 1' as ex P ecte d for unit total proba- 
bility. We then show that this chain distribution can be altered 
by further storing / subsequent qubits, resulting in the I th 
collective 1-spin down-flip, 



(9) 



fe=i 



with 



a 



min{;,fc-l} 

W(fc) = a 1) (fc) £ (-iyr Mr) 



r=0 



= a { Q \k) 2 F 1 (l-k,-iA;-t&n 2 9) cos 1 9 (10) 

where the more compact form 2^1(0, b;c;z) denotes the 
Gauss hypergeometric function JH. The detailed derivation is 
given in the Supplementary Material, together with the proof 
for decoding the stored qubits in the reverse order 01 ill . This 
also means the ferromagnetic chain can store the two-qubit 
states |t/i4/ 2 )' l4/it/ 2 )' which can be read out sequentially. 

We have also proved that encoding and decoding two \ 
qubits in the chain \F) C can be done in a similar fashion ll ill , 
with the (0, 0)*' 1 collective 2-spin down-flip amplitude (after 
encoding only 1-4 /1 4-/ 2 ) ) being 



a^L(fci,fc 2 ) = (■ 



(0,0) 



-isin(9) 2 cos fcl+fc2 " 1 



(2-(fc 2 -fci-2)tan 2 6 
(11) 

By linearity, the memory can thus store at least two qubits 
of arbitrary states, including any bipartite entanglement (and 
hence qubits of mixed states). A numerical check for N = 8 
static spins and two fully entangled flying qubits is shown in 
Fig. [3] demonstrating that this state can be retrieved with high 
fidelity and errors within the expected limits. 

To generalize the derivations, we need to show that the 
chain can store a spin state that includes an arbitrary number 
n of down-spins 



D 



(*",...,/„) _ X] a (;i!...,z„)( fcl '-"' fc ") %i,...,fc„) ( 12 ) 
fci<...<fe„ 



where the fc^s denote the spin-flip positions and denotes the 
number of \\) t encoded between the i th and (i + l) th flying 
|4) f qubits. This would mean that the chain can store a num- 
ber of flying qubits, each of which was originally either |f)j 
or |4) f and can be further retrieved by the aforementioned 
decoding mechanism. Arbitrary multi-particle superpositions 
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FIG. 3. Density state tomograms constructed for the corresponding, 
retrieved (two-qubit) quantum state from a memory chain \F) con- 
sisting of N = 8 static spins without decoherence. (a) and (b) repre- 
sent two separate cases of memory storage and retrieval of quantum 
information. The inputs, both entangled, are indicated in each case, 
together with the corresponding infidelity of the retrieved state. We 
chose a non-special value 9 = 1.1; for a tolerance e ~ 1CP 4 , this 
requires N > 6 to reliably store one qubit, as predicted by Eq.[7] 



and entanglement then follows by linearity. Unfortunately, 
this procedure becomes impractical for more than 2 down- 
spins due to the difficulty of keeping track of all the indices, 
and a rigorous proof for an arbitrary number of up-spin and 
down-spin qubits remains an open challenge. 

However, further support for the conjecture is provided by 
numerical simulations as shown in Fig.|4]for short chains with 
a few flying qubits of randomly generated states, and 8 chosen 
to give a high probability for a ferromagnetic tail after encod- 
ing. Analysis of the fidelity after decoding gives results which 
are consistent with the conjecture within the expected error 
bounds. 
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FIG. 4. Plot for the infidelity 1 — F of the retrieved qubit (with 
phase corrected), relative to the corresponding input pi, against the 
ordinal number i of the inputs. Here, four randomly-generated pure 
qubits {pi}i are sequentially encoded into the memory chain \F) C 
(N = 9), and then retrieved one by one in the reverse order without 
decoherence. We used 9 = 1, and hence for e ~ 10 -4 (or 10 -2 ) it 
requires N > 7 (or 4) to reliably store one qubit. The finite chain 
length restricts the number of qubits the memory can hold; within 
its capacity the last encoded (or the first retrieved) qubit always has 
the highest fidelity, since part of the quantum information stored for 
the earliest qubits may have been "pushed out" of the memory by the 
later ones during encoding. 



We conclude with a discussion of a potential architecture 
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FIG. 5. The static pseudo-spin structure, with a ground doublet plus 
a virtual, excited level: The XY exchange model between the flying 
photon and each QD can thus be achieved. 



for the memory. In essence, it is a slight extension of Jfle cav- 
ity quantum electrodynamics system suggested in @]. The 
qubits correspond to the polarization states of photons, which 
are transported reliably in a semiconductor photonic waveg- 
uide. An array of embedded quantum dots (QDs) act as the 
chain, where each multi-level QD exhibits an effective pseudo 
spin- 1/2 within its degenerate doublet ground state; an XY 
exchange model is achieved via two-photon Raman processes 
involving virtually populating an excited level (see Fig. [5j & 
A low operating temperature on the mK scale is required to 
initialise the memory, and also to have a long coherence time 
(~/xs) for each QD and hence the memory. With weak cou- 
pling between the photon and the QDs, the stored qubits are 
delocalised in the memory: the QDs store the qubits collec- 
tively. The QDs should also be well separated to have negligi- 
ble interactions with each other, say spaced by 10 to 100 nm. 
In that case, within one millimetre there can be as many as 
10 4 ~10 5 non-interacting QDs in the memory, where it only 
takes 10 ps for each write/read operation. In this optical struc- 
ture, 9 variations are also negligible, as are the variations in the 
kinetic energy of the photons. Therefore, this memory struc- 
ture is resilient to various imperfections. It should passively 
store qubits much more quickly, and have less photon loss, 
than previous holographic schemes since no active con- 
trol is needed during successive write/read operation rounds. 
Other systems such as optical lattices, where XY exchange 
interactions can be effectively achieved (and tuned), with an 
internal electronic degree of freedom of the cold atoms serv- 
ing as both the static spins and the flying qubits lioll . are also 
potential candidates for the memory realization. 
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READ-OUT OF THE FIRST QUBIT 

In the main text, we have shown that the quantum infor- 
mation held by the first flying qubit \4>i) ^ can be encoded 
into the memory chain via the proposed mechanism. The col- 
lective state of the chain after encoding is described by Eq. 4 
in the main text, while the flying qubit emerges as polarised 
It)/!- We further proposed that the state can be decoded 
back from the chain to a polarised flying qubit |t)/ injected 
with the same kinetic energy in the 'read' direction (see Fig. 1 
in the main text). Here, we shall derive the total state of the 
system after decoding (Eq. 5 in the main text), which provides 
further insight into the nature of the spin chain memory. 

Just before the read operation, the total state of the flying 
qubit and the chain is 



N- 



ai - ifii sin (9 ^ cosk ^ 6 $k ) It: F ) ■ 



fe=i 



(SI) 



The state |t, F) with amplitude a\ again remains the same 
during decoding, while each state |t, F), with amplitude 
— ij3\ sinflcos' 0-1 9, evolves as follows, as the flying spin 
passes each member of the chain: 

S k \tl F ) =^JV(->oo) ••• =>fc+l Sfc \t,F) 

=^ fe cosggfe |t,-F) -isin9 \i,F) 
=^>k-i cos 9 \t,F) - ism 9 x 

( cos 9 ||, F) - i sin 9 |t, F) ) 
=* fc i ... Vfc' = fc-2,...,2 
=>i C0S9S- \t,F) -ism9cos k - 1 9\i,F) 

k-1 

- sin 2 9 ^2 cos i_1 9 |t, F) if fc > 2 

or cos6SiVt,F)-ism6\l,F) if fe = 1 (S2) 

where we have followed the same notation as in the main text. 
Thus, by linearity, the total state after decoding becomes 

cxi |t, F) - ifix sin 6 ( cos 6 5f |t, F) — i sin 9\i,F)) 

OO 

ifix sinfl^cos^ 1 6>(cos(9 \t,F) — i sin 9 cos fc 1 9 |J,, F) 

k=2 

oo n—1 

tfx sin 9 cos' 1 - 1 0( sin 2 9 ^ cos^ 1 9 S~_ J |t, F) ) (S3) 

n=2 j=l 



where in the last line we have replaced the dummy variable k 
by n. By absorbing the bracketed terms in the first line into 
the summation of the second line, Eq. [S3]is equivalent to 



ti \t,F) -f3i sin 2 fl^cos 2 ^ 1 ) 9\i,F) J -iftsinfl 

k=l ' 



E cosk 6 $k 



k=l 



1 

oo n—1 

- 2 *EE 

n=2 j=l 



cos n+j - 2 9 S' 



\t,F), 



(S4) 



(i.e., Eq. 5 in the main text) where the first sum is simply a 
geometric series. For the double summation we want to focus 
on the coefficients for , and hence the sum over j can be 
replaced by a sum over k with k = n — j > 1, i.e., 



oo oo 



£ J2 ™ n+j - 2 Sn-, = E E cos 2 ^ 1 )-* 9 S k 

n=2 j=l k=l n=k+l 



fc=i • 



cos^ 9 



Si 



(S5) 



where we have swapped the double summations as the sum 
over n runs to infinity. Therefore, the second line in Eq. |S4] 
becomes zero due to complete cancellations for each k = 
1,2,... oo, and the total state after decoding is then simply 



(ai |t) - Pi \l)) f <g)\F) 



(S6) 



The complete cancellations only require that N — > oo, 
which is the same condition as for the encoding procedure. 
Thus in this limit, the state of the flying qubit with matching 
kinetic energy can be transferred to and back from the mem- 
ory in the way we have proposed in the main text. Note that 
this read operation is not a time reversal in general. 

If we employ the Heisenberg model, the state of the flying 
qubit can still be encoded into the chain, as explained in the 
main text. However, once read back, the total state of the sys- 
tem becomes (by tracking the extra phases e lk 8 in the above 
derivation), 



(ai |t)+#ll))/(g 
where 

oo 

p[ = -ft sin 2 6» ^ 

k=l 



l^) c + I>lt}/(g)^l F >c (S7) 
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7fc = -i/3ie iy sinf? ( e ike cos k 9 - 



sin f > e 

n=k+l 



e i(2n-fc)6 cos 2(« 



sin 6* cos 6* e 



1 - e 



2/7) 



1 - e 2ie cos 2 



(S7a) 



for all k e Z+ . The decoding procedure no longer works due 
to the extra phases, except when 6 = \ , Note that Eq. [S7] 
is consistent with the XY model once we identify each extra 
phase e' 10 with 1 accordingly. 

When the kinetic energy of the flying qubit during the read 
operation differs from that of the input for encoding (i.e., 
9dec 7^ 9 enc ), we numerically simulate how the error effects 
on the memory for various chain lengths TV in Fig. IS II where 
X = (&dec — 8 enc )/6 enc , Here, 9 enc is adjusted so that each 
chain with length N is just enough to collectively store the 



qubit (it is defined by setting cos 9 
Fidelity 



0.01). 







i 0.8 








/ 0.6 





-0.5 0.0 



0.5 



N=l 

N=5 

N=10 

N=50 

N=100 



FIG. SI. Plots of the fidelity of the retrieved qubit, relative to the 
input |+) . , against the fractional difference \ between the coupling 
strengths during encoding (9 enc ) and decoding (#dec), for various 
chain lengths. For each N, the whole chain collectively stores the 
qubit, i.e., N 3 = N. 



From Fig. IS 1 1 we see that after the decoding round, the re- 
trieved qubit is of a high fidelity (> 99%) with respect to the 
original input, for small mismatches (~40%) in the encod- 
ing and decoding 9 values. Moreover, as the number N of 
static spins increases, the memory's tolerance to such errors 
improves. 



COLLECTIVE 1-SPIN DOWN-FLIP DISTRIBUTIONS 



In the main text, we introduced the I th collective 1-spin 



down-flip distribution D 



= H k ^\k)S^ inEqs. 8-10; 
(A;) | 2 corresponds to the probability of the k th static spin 
being \l) Sk in the respective distribution. 

Arising from the chain state after encoding the first flying 
\i) qubit (see Eq. 3 in the main text), the th 1-spin down- 
flip distribution has \k) = —ism9cos k ~ 1 9. The total 
probability is J^'kLi \% Wl 2 = 1' while the mean position 



/Xq of the down-flip and the associated standard deviation 
cr^ are (see Appendix I) 

oo 

^fc |4 1} (fc)| 2 = esc 2 6» 



fc=l 



^fc 2 |4 1} (fc)| 2 - (4 1} ) 2 -cosflcsc 2 ^. (S8) 

\ k=l 



The I th collective 1-spin down-flip distribution in the mem- 
ory results from further encoding I subsequent flying |f) . 
qubits. We now inductively derive the analytical expression 
for a\ (k) given by Eq. 10 in the main text. 

First, we find that after encoding one flying qubit into 

a chain with state S k \F) c V finite k < N, the total state 
becomes 



\t) f (g){co S 6S^-sm 2 6 ^ cos"-^ 1 8 S~)\F) C 

n=k+l 

(S9) 

which is derived step-by-step as before. Thus, by linearity 
(apply Eq. |S91 l, encoding one flying qubit into the chain 

with distribution Dq = J2 k a^\k)S k gives rise to the first 
1-spin down-flip distribution 

oo oo 

D[ 1] = ^ 4 1) (fc)( cos 6*^7 -sin 2 6» ^ cos™^ 1 9 S~) . 



k=l 



n=k+l 



(S10) 

Again, we want to focus on the coefficients of S k terms for 

the double sum. Note that \k)/ cos fc_1 9 = —isin9, in- 
dependent of k. Thus, the double sum in Eq. lSlOl is 



oo oo 



sin2e E E 4\k) C os n - k - l 6S- 



oo oo 



= -tan 2 (9^] (-isin6)cos n eS- 

k=l n=k+l 

oo oo 

= -tan 2 (9Hsin(9)^ ^ cos fc 8 S k 

oo 

= - tan 2 9{-ism9) ^T(k - 1) cos fe 8 /S~ 
fe=i 

■oo 

= -tan 2 9 E( fc ~ !) a Q ) ( k ) cos 6 % ( S11 ) 



fc=i 



where we have interchanged the labelling of dummy variables 
n and k in the third line, and evaluated the sum of geometric 
series in the fourth. Substituting back into Eq. IS 101 (and the 

more general definition of D\ f 1 '), we have 

a^\k) = 4 1} (k) cos 9 (1 - (k - l)tan 2 <9). (S12) 

We can then apply the same procedure (with Eq. |S9b to 
encoding one flying If)/ qubit into d[ 1} = J2k a i\ k )Sk 



3 



to find, again by linearity, amplitudes for the second 1-spin 
down-flip distribution 



4 X) (fe) = a^\k) cos 2 9 x 



(l-2(fc-l)tan 2 fl + {k 1} o (fc 2) — 4 ' 



tan 



(S13) 



where (as also illustrated in Eq. IS1 II ) we have evaluated the 
following weighted sums of geometric series (by first rela- 
belling the dummy variables to focus on SZ), 



oo oo 



oo oo 



k=l n=k+l 



n— 1 k— n+1 

OO ✓ fc— 1 \ 

oo 

= E 



fc=2 N n=l 



^y-^cos^r. (S14) 



fc=3 



Note that the general term of the bracketed series are obtained 
from evaluation of the weighted geometric series in the previ- 
ous step (evaluating the lower order distribution). 

Now, with the key inductive steps to relabel the dummy 
variables (to focus on SJ^), and to evaluate the following series 
(as done above and in Eq. lSl U , 

fe-l l' -l l' , w' + i I'+l 



n=l' + l m— j — 1 

we can inductively derive 



cause no spin flips at all, or can move a spin down from a site 
nearer the front of the chain to one further along it. The th 
term results from the situation in which the down spin is ini- 
tially localized on spin state k and where all I subsequent |f)^ 
qubits move along the chain without executing further flips. 
Each has contributed a factor of cos 9 due to the exchange in- 
teraction and thus an overall factor of cos' 6 is present in addi- 
tion to a^\k). In general, the r th term (r > 0) occurs when 
r movements of the initial spin down position occur before 
that spin down reaches its final position k. There are ( ) ways 
of choosing the r qubits which cause the flips from the I total, 
and ( _1 ) ways of choosing which r of the (fc— 1) static spins 
which precede the k th will hold the spin down at some point 
before the spin down finally occurs at site k. The other terms 
in the summation come from the fact that each double spin flip 
(or movement of the spin down location) gives rise to a factor 
of (— isin6*) 2 , while losing a factor of cos 6*; in addition, the 
other (I — r) [f}* qubits passed the down-flipped spin without 
exchanging, and each contributed a factor of cos 9. Combin- 
ing these coefficients gives rise to Eq. IS 161 

Decoding this more general memory state can be achieved 
by injecting successive |j) , spins in the decoding direction. 
After the first such spin passes, the new memory state down- 
spin amplitude for the k th site, (fc), results from two pos- 
sible scenarios: Either this |f) , passed the k th site which was 
already in the down state, without exchange, or it transported 
the (k + s) th down spin to the k th position. Taking into ac- 
count of the factors contributed, we have 



(S 15) af*l (k) = 0; (1) (k) cos e+(-i sin 9) 2 E a ; W + a) cos 3 " 1 1 



a^\k) = a\ L, {k) cos' 9 x 



a«(/c). 



(S17) 



min{;,fe-l} . , — 1 . r 

E i-l) r hU{l-n)hUQ°-3)ta* ar 



r— m— j — 1 



To establish this last equivalence, we multiply both sides of 
the following identity (see Appendix II) 



min{/,fc— 1} 



a£\k) cos 1 9 E (-!) r 



r=0 



k - 1\ (I 



4 1} (k) cos 1 9 2 F 1 (l-k,-l;l;- tan 



r I \r 

2 



2-Fi(a, b;l;z) z 



tan 



1 - z 



1 - z 



E 

s=l 



2 F 1 (a - s,b; l;z) 
(1 - z)> 



2 F 1 {a,b+l;l;z) 
(SI 8) 



(S16) by a^\k) cos 1 1 9, and substitute for a = 1 — k,b 



Note that Eq.[S9]has two terms, the first (and lower order) of 
which inductively adds to the lower order terms (in tan 2 9) 
in Eq. IS161 to give rise to the bracketed product coefficient 
concerning I. □ 

In the more compact form of Eq. IS 1 61 2^1(0, b; c; z) 
denotes the Gauss hypergeometric function and 
ag (k) cos 1 9 renders possible divergence of tan 2r 9 conver- 
gent in Eq. IS 1 61 (analytic continuation is assumed implicitly 
here). The expression for the I th collective 1-spin down-flip 
distribution can also be obtained through a combinatorial ar- 
gument, as follows. For any fixed down-flip position k, and 
assuming k > (I + 1), the amplitude is a sum of / + 1 terms 
, each of which corresponds to a different origin for the . 
In general a spin-up qubit passing along the chain can either 



and z = — tan 2 9; we then obtain a^\{k) = a[]} 1 (k) from 
Eq. IS 171 Here, I £ N can be arbitrary. This means that the 
1-spin down-flip distributions can be manipulated in both di- 
rections, essential for the chain to act as a memory. 

Having established this important feature, we now show 
that each distribution corresponds to a unique storage mode 
and the modes are independent, i.e., expressed as the follow- 
ing (discrete) orthonormal condition 



E4 1} < 

k=l 



(k) aY\k) 



(S19) 



This ensures the unit total probability for each distribution. To 
establish this discrete orthonormal condition Eq. IS 191 we first 
introduce the normalized Meixner polynomials fljj] (with j, x 



4 



integers) 

M'Ax\ (i,v): = 2F1 -x; fj,; 1 - - 



Population 
0.4 



= %Fi y-x, -j; - -j (S20) 

where the equivalence comes from the symmetry of the hy- 
pergeometric function in its first two arguments. Note that the 
different normalization is present since we are only summing 
over x € N (instead of Z). The orthonormality condition for 
the Meixner polynomials states 

00 

^ Mj(x; fx, u)M'j, (x; /i, v)uj{x; /i, v) = 5jj> (S21) 
x=a 

where the discrete weight ui(x; /i, v) = (1 — i/)^-^ 2 -!^ Ill- 
Setting x = k — 1, j = I, fi = l,i> = cos 2 9, and substituting 
Eq.|S20]into Eq. lS2T1 we have 



]T 2 Fi(l-fc, -Z;l;-tan 2 <9) 2 -Fi(l - fc, 1; - tan 2 0) 
x sin 2 cos 2(k+l - 1] = S w (S22) 



fc=i 



which is exactly the desired condition Eq. IS 1 91 Note that one 
special case of the Meixner polynomials are the Krawtchouk 
polynomials which have recently been applied to works 
involving quantum state transfer of a single spin excitation 
within certain linear, interacting spin chains 12J]. 




FIG. S2. Plots with a non-special coupling strength 9 = 0.4 for (a) 



the mean positions fij 1 ' of the down-flip against the number I of sub- 
sequent flying |f) . qubits encoded into the memory; (b) the corre- 
sponding standard deviations <xf 1 ' against /. The linear relationships 
are also observed for other 6 values. 



The independence of these unique storage modes are best 
illustrated by considering their mean positions fij 1 ' of the 

down-flip and the corresponding standard deviations <r. , 
which can be calculated in a similar way to how we found 
Eq. |S8] We find that for any given coupling strength 9, both 
fjq and increase linearly with /; an example is demon- 
strated in Fig.|S2] with a non-special value 9 = 0.4, say. 

Finally, we numerically simulate (see Appendix III) the 
down-flip distributions, by encoding a flying \l) ^ qubit fol- 
lowed by I subsequent |f)j into the initially ferromagnetic 




FIG. S3. Plots of the down-flip distribution \a^ (k)\ 2 in a chain with 
N = 9, and 6 = 1.2: The discrete plots are from our simulations 



("blue circle" for la^'COl^ ^ " re d square" for |t£ '(k)f), while 



'(i), 



the curve corresponds to the analytical solution Eq. lS16l (onlv values 
for integer k are relevant). 



chain. An example is shown in Fig. [S3] and shows agreement 
with our analytical solutions for |tij (k)\ 2 . 



THE (0, 0) th 2-SPIN DOWN-FLIP DISTRIBUTION 

To generalize the results, we need to show that the chain 
can store a spin state that includes an arbitrary number n of 
down-spins 



D 



(n) 

(h,-M 



~ a {l u ...,l n )(klT--ik n ) S(ki,...,k n ) 



fci<...<fc„ 



(S23) 

where the fcjS denote the spin-flip positions and Zj denotes the 
number of |f)^ encoded between the i th and (i + l) th flying 
qubits. This would mean that the chain can store the 
information from a number of flying qubits, each of which 
was originally either |f)j or \\) * and can be further retrieved 
by the aforementioned decoding mechanism. By linearity, it 
could also store any superposition, which would confirm its 
status as a true quantum memory. We, however, are unable to 
prove the general case due to the increasing complexity of the 
analytical solution (with large numbers of parameters Z;S and 
kjS). 

Here, we derive the analytical expression for the (0, 0) th 
collective 2-spin down-flip distribution 

5 (o. ) o)= E a (o!o)( fc i^)5- lifc2) (S24) 

ki<k 2 

after having encoded only two in the chain. Going 
through step-by-step we find that the total state of the system, 
after encoding one flying \l) ^ qubit into a chain with state 

$k \F)c> becomes 



|t) / 0(-iBinfl) J2 co S k '~ 2 6S^S^ \F) C iffc = L 

fe'=2 

(S25) 



5 



or V finite 2 < k < N, 
lt>/0(-isiii< 



fc-i 



^2 cos k '~ 2 e s-s-, + ^ cos fe ' e s~s k 

■k'=k+l k' = l 

k—1 oo 



tan 2 9 y~] cos 

k' = l k"=k+l 



k'+k"-k 



^k' Sk" 



\F)c 



(S26) 



Note that by applying a combinatorial argument as illustrated 
before, one can obtain the same results. 

Thus by linearity, after encoding the \l) ^ into Dq \F) c , 
we have 



D 



(2) 
(0,0) 



(— i sinl 



r oo oo 

E E 

fc=i k'=k+i 

OO / k — 1 

e(e 

fe— 1 oo 



COS 



fc+fc'-l 



tan 2 9 J2 E cosfc ' +fc "~ 16 '^^ 



fc' = l k"=k+l 



(S27) 



By manipulating the double and triple summations as before, 
we find that 



four randomly-generated, pure qubits, with the same coupling 
strength (9 l = 6 V i). Here, we consider the effects of the 
different kinetic energy that each qubit may possess, and the 
memory's tolerance to the resulting coupling strength varia- 
tions. 
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FIG. S4. Plot for the infidelity 1 — F of the retrieved qubit (with 
phase corrected), relative to the corresponding input pi, against the 
ordinal number i of the inputs. Here, four randomly-generated pure 
qubits {pi}4 (same as in the main text) are sequentially encoded into 
the memory chain \F) (N = 9), and then retrieved one by one 
in the reverse order without decoherence. 6 l denotes the coupling 
strength between the i th qubit pi and each static spin during both en- 
coding and decoding rounds; site-to-site variations are ignored here, 
and considered later. Note that for (b), ten independent runs of simu- 
lations are performed, and for each qubit the average is taken for the 
phase-corrected fidelity (to discount the imperfect randomness in 0' 
for small number of rounds n and hence possible outliers). 



D 



(2) 
(0,0) 



-I sine 



E E 

ki — 1 ki— fei+1 



COS 



fcl+fc2-l 



,fel+fe 2 -3 , 



- (fc 2 - fa - 1) tan^ 9 cos fcl+fc2 - 1 9 ) S£ S£ . 

(S28) 



Therefore, the (0, 0) collective 2-spin amplitude is 



j (2) 

(0,0) 



(fa, fa) = (—i sin( 



fci+fe 2 -i 



1 (2-(fa-fa-2) tan 2 9), 
(S29) 

as quoted in the main text. We have the total probability 

Sfc 1 <fe 2 l a (o o)(^ 1 ' fa)\ 2 = ^ ne Qubits can also be recov- 
ered sequentially, as can be shown by going through step-by- 
step or using a combinatorial argument. 

We have thus shown that the memory can encode and de- 
code It/st/i), Whih)> \ih t/i>> and \lhlhY> b y linearity, 
the memory can store two qubits of arbitrary state, entangled 
or not (and hence qubits of mixed states). We conjecture that 
the memory chain can store multiple flying qubits of arbitrary 
states using the described mechanism (see an example of nu- 
merical simulations of larger memory states from Fig. 4 in the 
main text supporting this conjecture). 




(a) Pi 




(b) PT 



FIG. S5. Density state tomograms ("red" for real components and 
"yellow" for imaginary) constructed for (a) the input pi generated 
above (i.e., from the main text); (b) the retrieved qubit p\ et , by 
encoding only p\ into the the memory chain \F) (with N = 9) 
and then read back (without correcting the phase). We have random 
9k £ (0.9, 1.1) between the qubit and the k th static spin. The infi- 
delity of the retrieved qubit, if phase corrected, is 7.4 x 10~ 6 , which 
is again averaged over the results from ten independent runs of sim- 
ulations. For comparison, the corresponding infidelity for 6k — 1 Vfe 
is 6.1 x 10" 6 . 



VARIATIONS IN COUPLING STRENGTHS 

In the main text, we have shown an example of numeri- 
cal simulations involving the memory storage and retrieval of 



To separate the effects of finite chain length, we employ 
the same qubits {pi}4 generated in the main text and simu- 
late their storage and retrieval from the memory \F) with 
the same length N = 9, but with random 9 % 6 (0.9, 1.1) for 
each i. Each 9i matches for the corresponding decoding and 



6 



encoding operations of pi (effects of such mismatches were 
considered in Fig. IS lb ; site-to-site variations are also ignored, 
but considered below. Fig. [S4]thus illustrates the effects of 
small round-to-round variations (~ ±10%) in the coupling 
strengths, and shows that the memory is robust towards such 
imperfections as long as the memory has sufficient capacity to 
store each qubit reliably (i.e., has a large enough N). 

To consider site-to-site variations in 6k, we encode and de- 
code only pi into the chain \F) with N = 9, with variable 
coupling strengths 6k G (0.9, 1.1) between the qubit and k th 
static spin. Note that again 6k matches for the decoding and 
encoding operations for the same position k. In this case, we 
found the infidelity (averaged over ten independent runs) be- 
tween p\ et and pi to be ~ 10~ 6 , similar to that in the case 
of 6k = 1 Vfc (see Fig. IS5I ). Thus, the memory can also tol- 
erate small site-to-site variations (~ ±10%) in the coupling 
strengths 6k- 



POTENTIAL ARCHITECTURES 



APPENDIX I 

Calculations of the mean and standard deviation in Eq. 
involves evaluating series of the following forms 



OO 7 OO 7 / 

E*~-£(IV)-s(rf 



fe=l 



fc=l 



dy V 1 - y 



E* 

fc=i 



X -E 

fe=l 

d?_ 

~ dy 2 
dy 2 



(k + l)ky 



k-l 



oo 

E 

k=l 



ky 



fe-i 



fe=i 



1 



i 



yf 
-y 



{i-yf 



(i - yf 



(A2) 



We then set y = cos 2 6 for the two final formulae, and substi- 
tute back into Eq.[S8]for evalutions. 



APPENDIX II 

In order to show that a'^\ (Eq. IS17I ) is the same as a^ v 
we applied the identity stated in Eq. IS18I 



2 F 1 (a,b;l; z) 



1 



1-z 



^ 2 F 1 (a- s,b;l;z) _ 

E (T^zy = 2^1(0,6+1;!; z). 



In the main text, we have discussed the semiconductor pho- 
tonic waveguide system JH], suitable as a candidate for our 
memory system. We have also suggested the more specula- 
tive optical lattice, which provides the advantage of tunable 
exchange couplings by controlling the laser fields fl. A fur- 
ther alternative could be a solid-state system of electrons mov- 
ing in a gated semiconducting quantum wire. The memory is 
a chain of QDs, each accommodating an electron spin. A fly- 
ing qubit is manifested in the spin of a moving electron. A 
drawback of this system is that the flying-static qubit-spin in- 
teraction is of the Heisenberg type |5], and thus one needs to 
work with 6 ~ ^. The storage of quantum information is 
between the same species, and is localised in the memory. 

It remains a speculation (and requires further investigation 
to see) whether our proposed memory system could be ex- 
tended to perform full quantum computation iQ-Ht] . 



To prove this, we work with the Euler's integral representation 
for the hypergeometric function 

- 1 t b - 1 (l-t)- b 



2 Fi(a,b; l;z) 



1 



-dt (A3) 



B(b,l-b)J (l-zt) a 

for \z\ < 1, and apply analytic continuation to other values of 
z. Here, B(p,q) = T(p)T(q)/T(p + q) is the beta function, 
where T{p + 1) = pT(p) is the gamma function 01. We see 
that 



B(b, l-b) = 



r(b)m-b) |r(6±l)(-6)r(-b) 

r(i) 



r(i) 

= -B(b+l,-b). (A4) 

Thus, by substituting both Eqs.[A3]and[A4]back into Eq. lS18l 
(and multiplying it both sides by (1 — z)B(b, 1 — b)), our task 
reduces to proving that 

t b - 1 (l-t)- b J ^ z f 1 t b -\l-t)~ b 



/' 

Jo 



j=i (1 - *) s Jo (1 - zt) a 



-dt 



r 1 t b (i-t)- b - 1 , 



(A5) 



jo {i-zty 

where the extra minus sign on the right hand side comes from 
Eq. IA4I Evaluating the summation in Eq. IA5I we have 

v ' -dt 



-dt 



£j (i - z y j (i - zty 



1-Zt 



si j.b-1 



t b ~ L (l-t)- 
(1 - Zt) a 



(Al) 



l x j i-zt \ t b -\i-ty b 
~ I ?{/(t-i)J {i -zty dt - 



(A6) 



7 



Hence, the left hand side of Eq. lA5l is equivalent to 

t (i 1 - zt Y~ l{1 ~ t) ~ h dt 

J a V i-*y {i-zty 

= -(1-^) / ,/ ' dt, (A7) 

which is exactly the right hand side of Eq. IA5I □ 



APPENDIX III 

The time evolution operator for the interaction between a 
flying qubit and a static spin can be expressed as in Eq. 2 in 
the main text. However, each flying qubit interacts with the 
static spins sequentially; in order to numerically simulate the 
evolution for the state of the whole system \ifjf, <& c )> we a Pply 
the write and read operators, 

U wr ite = UfN---Uf2Ufi 
Uread = U fl Uf 2 ...U fN (A8) 



where l m is the m x m identity matrix, C m = cos Ok l m 
while S m = — ismOk l m ; both m and blank fields cor- 
respond to blocks of Os. The dots represent the alternating 
repetitive patterns (along the diagonal in the upper left and 
lower right blocks these are t m , C rn , l m , C m along the 
diagonals in the lower left and upper right blocks these are 
S m , 0m> S m , (D rn ...), and there are an equal number of blocks 
of l m 's, C m 's and S m 's. This is readily obtained by con- 
sidering the exchange processes between the flying qubit and 
the k th static spin, regardless of the state of other spins (the 
alternating repetitive pattern due to spins s\ to Sk-i and the 
size of each block due to Sk+i to s«). Note that each Ufk is 
symmetric, and centrosymmetric (symmetric about the centre 
point and a property that is closed under matrix multiplica- 
tions) H. 

When the whole system is restricted to at most one excita- 
tion (as we did in the main text for considering decoherence), 
the state space consisting of more than one excitations is never 
accessed and can thus be ignored for the purpose of simula- 
tion. In this way, the dimension (N + 2) x (N + 2) of the 
density matrix (with only zero or one excitation) grows lin- 
early with the number N of chain spins; simulations of very 
long chains can then be done conveniently. 



i th static spin 



where each 2 N+1 x 2 N+1 time evolution matrix U fi describes 
the interaction between the flying qubit and the 
under the standard basis {/, si, s 2 , sat}): 



u Ik = 



) 2 N- 



c 



2 N- 



J 2 N- 



1 2 

(A9) 
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